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Gravitational radiation reaction in compact binary systems: 
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The radiation reaction in compact spinning binaries on eccentric orbits due to the quadrupole- 
monopole interaction is studied. This contribution is of second post-Newtonian order. As result 
of the precession of spins the magnitude L of the orbital angular momentum is not conserved. 
Therefore a proper characterization of the perturbed radial motion is provided by the energy E 
and angular average L. As powerful computing tools, the generalized true and eccentric anomaly 
parametrizations are introduced. Then the secular losses in energy and magnitude of orbital angular 
momentum together with the secular evolution of the relative orientations of the orbital angular 
momentum and spins are found for eccentric orbits by use of the residue theorem. The circular orbit 
limit of the energy loss agrees with Poisson’s earlier result. 


I. INTRODUCTION 

Compact binary systems composed of neutron stars and/or black holes are among the best candidates of sources 
emitting gravitational radiation in the frequency range of the Earth-based interferometric detectors such as the Laser 
Interferometric Gravitational Wave Observatory (LIGO) 0, VIRGO (2|, GEO |jj, and TAMA Q, and also for the 
envisaged Laser Interferometer Space Antenna (LISA) . Before the final coalescence of the compact binary due 
to a cumulative loss in energy and angular momentum there is a regime where a post-Newtonian description of the 
motion and of the gravitational radiation is well suited, provided it is taken with sufficient accuracy. 

As up to the second post-Newtonian (PN) order the energy and total angular momentum of the system are conserved, 
and the leading contribution of gravitational radiation at 2.5 PN is generally agreed to decrease the eccentricity of 
orbits faster than their radii |t|, it is customary to compute radiation templates and the radiation reaction with the 
assumption of circular orbits. Another argument for such a simplifying assumption is the reduction of the number of 
fitting parameters in the matched filtering. However as described in Refs. (Js| and |J, there are various astrophysical 
situations in galactic nuclei where the eccentricity of the orbits is still significant. The error one makes in searching 
for waves emitted by binaries on circular orbits when the waves come from binaries on eccentric orbits was estimated 
to be substantial in JT(^]. Therefore it is desirable to have a generic treatment, also valid for eccentric orbits. Such a 
treatment, valid up to the second post-Newtonian order was provided by Gopakumar and Iyer ]lT| . 

But the situation is even more complicated. At 1.5 PN order other features of the binary appear. These are spin- 
orbit type contributions, which further increase the number of parameters to be fitted. At 2 PN both spin-spin and 
quadrupole-monopole effects appear. The latter are the consequences of considering each component of the binary as 
a monopole in the quadrupolar field of the other. The nonradiative dynamics of spinning binaries with emphasis on 
both the spins and the quadrupole-monopole interaction was investigated long ago in a series of papers by Barker and 
O’Connell (see, for example, iQ). Radiative aspects concerning spins were studied in Refs. |Ti|- However the 
quadrupole-monopole type radiative contribution has received less attention until now, excepting the circular orbit 
limit, discussed in detail by Poisson po| . In this paper we extend this treatment to eccentric orbits. 

For this purpose in Sect. II we derive two parametrizations of the orbit, which generalize the Keplerian true and 
eccentric anomaly parametrizations for the case of the quadrupole-monopole perturbation. These parametrizations 
have the advantageous property that the required integral expressions can be easily computed by use of the residue 
theorem and, further, in the majority of cases the only pole is at the origin ©• Our treatment follows closely and 
relies on the results of Refs. Q and |li|. Here we also find the peculiar feature encountered previously solely in the 
discussion of the spin-spin contributions, that the magnitude of the orbital angular momentum is not a conserved 
quantity. We bypass this difficulty by introducing its angular average L and characterizing the perturbed orbit by 
E and L. It is equally possible to characterize the motion by E and the time average (L); however the resulting 
expressions would be more cumbersome. For completeness, the relation between ( L ) and L is given. At the end of 
the section we find the same functional expression for the period of the radial motion as in the unperturbed case, 
but with the energy characterizing the perturbed motion. This is a consequence of the specific functional form of the 
perturbing term in the radial equation |HJ. 

Having developed the required toolchest, the computation of the quadrupole-monopole contributions to the losses 
of energy and magnitude of orbital angular momentum together with the evolution of the angle variables, which 
characterize the relative orientation of the spin vectors and orbital angular momentum vector are straightforward. 
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These are given in Sect. Ill as the main results of the paper. We will show in the concluding remarks that in the 
circular limit the energy loss reduces to the expression given previously by Poisson. 

The velocity of light c and the gravitational constant G are kept in all expressions. 


II. THE RADIAL MOTION 


The Keplerian motion of the binary perturbed by the monopole-quadrupole interaction between its components is 
properly characterized by the Lagrangian 1201, p2||: 


C — Cn + Cqm , 

/xv 2 Gmp 
L-n = —x—I- 


Cqm — 


Gpm 3 
2 r 5 


z r 

^Pi 3(s ; -r) 


( 2 . 1 ) 


where Cn represents the Newtonian part of the Lagrangian (with m the total mass and /t the reduced mass). The 
summation is taken over the two bodies in the quadrupole-monopole interaction term Cqm > which depends both on 
the direction Si of the spin of the i th body Si and on pi, defined as 


Pi = 


Qi 

rriim 2 


( 2 . 2 ) 


Here Qi is the quadrupole-moment scalar | 
In the Newtonian limit computation gives 


of the r 1 axially symmetric binary component with symmetry axis Si. 


O = 0' — 0 = _ — (— — 1 
Ql 1 1 fli 10' 


(2.3) 


for each axisymmetric body characterized by the principal moments of inertia (0,;, 0j, O') and angular velocity 
Cli = Si/Q'i. 

Other contributions, like PN, 2PN [;fl|, spin-orbit (SO) |l7| and spin-spin (SS) Q, |p| appearing at 1.5 and second 
post-Newtonian orders, respectively will add to the dynamics. However these other post-Newtonian terms will not 
interfere with the quadrupole-monopole interaction terms. Quadrupole-monopole terms are singled out from among 
the rest of the second post-Newtonian order terms by the distinguishing parameters pi in the same way spin-spin 
interaction terms can be recognized as being quadratic in the spin magnitudes. Therefore up to the second post- 
Newtonian order each contribution can be computed independently, including the quadrupole-monopole one. In fact 
the first interference terms appear in t he 2 .5 order, as the first PN correction of the SO-terms [p3| . 

The acceleration resulting from Eq. (2.1) is 


a = & N + a q M 
Gmr 

&-N q 5 


a QM = 


3 Gm 3 
2 r 7 

3Gm 3 


£ r 2 

^Pi 5(s ; -r) 
i= 1 L 
2 

(§i • r) Si . 


(2.4) 


i =1 


By performing the required changes in notations the acceleration a qm coincides with Eqs. (55)-(56) of Ref. @- 
The Lagrangian governing the quadrupole-monopole interaction is not time dependent; therefore the total energy 
E of the system is conserved. We further remark that, as it happens for spin-spin interaction, here the Lagrangian 
Cqm has no dependence on velocities either. Therefore p = /tr and two conclusions emerge: 

(a.) The energy is given by 


E = En — Cqm , 


(2.5) 
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(b.) There is no quadrupole-monopole contribution to the orbital angular momentum 


L = r x p = L w . (2.6) 

The functional form of the Newtonian energy and orbital angular momentum vector flg| allows one to relate them to 
v 2 and r 2 : 


o 2 Em 2 Gm 

' = -+ - , 

fi r 

2 2 Ejy 2 Gm L 2 n 

/i r fj?r 2 


(2.7) 

( 2 . 8 ) 


The cross product of Eq. (2.4) with pr gives the evolution of the orbital angular momentum: 

3Gpm 3 


L =- 




i= 1 


r r x Si 


(2.9) 


Due to the quadrupole-monopole interaction the spins undergo a precessional motion, computed a long time ago by 
Barker and O’Connell [Eqs. (39) and (43) of 0 referring to bodies with arbitrary, but constant mass, spin and 
quadrupole moment]: 


Si = - 


3 Gpm c 


■Pi(Si ■ r) (r x Si) . 


It is then straightforward to check the conservation of the total angular momentum 

J = L + Sx + S 2 . 


( 2 . 10 ) 


( 2 . 11 ) 


Multiplying Eq. ( |2.9| ) by the direction of the orbital angular momentum L the magnitude of the angular momentum 
is found to evolve. The nonconservation of the magnitude of the orbital angular momentum is a feature already 
encountered when studying spin-spin dynamics. By performing the transformations enlisted in Ref fig)] the detailed 
form of the evolution of L emerges 


L 


3 Gpm 3 
2 r 3 


2 

y Pi sin 2 m sin 2 (ip - ipi) ■ 
2 = 1 


( 2 . 12 ) 


Here and hereafter we employ the following angles, k,; = cos _ 1 (Si-L) and 7 = cos _1 (Si - 82 ) characterize the 
relative orientation of the angular momenta. The angles ip and ipi are subtended by the intersection line of the planes 
perpendicular to J and L with the position r and the projections of the spins in the plane of the orbit, respectively 
[see Fig. 1 in Ref. Q. We also introduce Si = 2 (ipo — ipi)\. This description allows for computing the value of L 
in terms of the Newtonian true anomaly parameter % = ip — ipo and a properly defined angular average L of L in a 
similar way as was done in the spin-spin case: 


L( X ) = L + SL 


5L 


^4 jxT ^ Pi sin 2 Ki {2A cos(x+Si) 
2—1 

+ (3Gfj,m+2Acosx) cos(2y+i5,:)} . 


(2.13) 


where 


A = ^G 2 m 2 p 2 + 1/2 (2-14) 

is the magnitude of the Laplace-Runge-Lenz vector for a Keplerian motion characterized by E and L. 

By following the method of Ref. |L 8 | the quadrupole-monopole interaction part Eqm of the energy can also be given 
in terms of x an d r '■ 
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Gfim 

E QM (r, X ) = -2^ 


Pi 


1 — 3sin 2 Ki cos 2 IX+ 77 


Then the expressions (2.7) and (T 8 ) take the form: 


y 2 = 2 [E - E QM {r,x)} 


= 2 [E - E QM {r,x)\ 


2 Gm 
r 

2 Gm 


L(X ) 2 

p 2 r 2 


(2.15) 

(2.16) 
(2.17) 


with L(x) and Eqm{t , x) given by Eqs. fl2.13| ) and (2.15). As the parameter x appears only in second post-Newtonian 
terms, it can be replaced by any generalized true anomaly parameter having the proper Newtonian limit. In what 
follows we will seek for such a suitable parametrization r = r (x) of the radial component of the perturbe d mot ion. 

For this purpose we compute the turning points of the orbit, defined as roots of the radial equation ( 2.17 ) with 
f 2 = 0 for x = 0, 7 r. We find: 


Gmp ± A Gp 2 m 3 


-2 E 


4AL 2 


> 


i=l 


p \p = <Xq(A -f Gmp) + /3q(4A 3 Gmp) , 

= 2 [l — 3sin 2 Kj cos 2 (t/>o —t^i)] , 

/?o = sin 2 cos ■ 

Then, following Ref. (l 8 | we derive the generalized true anomaly parametrization 

L 2 Gp 2 m 3 


r = 


p(Gmp + Acosx) 4AL 2 (Gmp + A cos x) : 


Em 1 


(2.18) 


(2.19) 


with 


A* = A[A 2 (a(,+4^) + (GTO/r) 2 (3a( ) + 10^)] 


+ G?n/r[A 2 (3ao + ll/3o) + (GTO/x) 2 (Q:o-l-3 / 9Q)] cosx • 

In a similar way we derive the generalized eccentric anomaly parametrization 


r = 


Gmp — A cos £ Gp 2 m 3 


-2E 


4 AL 2 


E 


p»- 


with 


= A(a o + 4/3q) + Gmp(a l Q + 3/3q) cos£ . 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


Note that p \ F , A®, and 2® have the same structure as in the spin-spin case; however, the angular expressions a l 0 and 
(3q are different from ao and /3q introduced in Ref. fll]. 


In the perturbative terms we will need the Keplerian relations 


r = -j- sin x , 


i> = -4 • 

pr - 


(2.23) 


The two parametrizations are suitable for averaging radial expressions by use of the residue theorem, as described 
in detail in [Ell. The time average of L(x) is 


(L) — L + 


Gm 3 p 2 F 1 
4 A 2 L 3 F 2 


2 

pi sin 2 Kj cos Si 

i=1 


(2.24) 


with the coefficients T \ j2 given by 
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Fi = 2L(-2/x£) 1/2 [A 6 - 15G 2 m 2 v 2 A 4 + 32G 4 m 4 /j, 4 A 2 - 16G 6 ro 6 /z 6 ] 

+Gmn 2 [-llA 6 + 58 G 2 m 2 f j 2 A 4 - 80 G 4 m 4 n 4 A 2 + 32G 6 to 6 /z 6 ] 

F 2 = 4GmL(—2fj,E) 1 / 2 [A 2 - 2 G 2 toV] + A 4 - 8 G 2 m 2 ^i 2 A 2 + 8G 4 mV . (2.25) 


We remark that the coefficients are the same as the correspo nding coefficients in the spin-spin case. As L (%) 
given in terms of ( L ) has a more cumbersome expression than Eq. ( [2. 13 ), we will give all forthcoming expressions in 
terms of L. 

For the period once again we find a Keplerian expression, but with E characterizing the perturbed dynamics: 


T = 2TtGm 



(2.26) 


III. LEADING ORDER QUADRUPOLE-MONOPOLE CONTRIBUTION TO THE EVOLUTIONS OF 
THE DYNAMICAL VARIABLES UNDER RADIATION REACTION 


A. Energy loss 

To leading order the instantaneous loss of energy under radiation reaction is given by Einstein’s quadrupole formula 

(3.1) 


^_ = _ j{3)jl j{3)jl 

dt 5 c 5 


where is the completely antisymmetric Levi-Civita symbol, the numbers in parantheses denote a corresponding 
order time derivative and the system's symmetric trace-free (STF) mass quadrupole moment I^ 1 to leading order is 
given by 


In = M {x 3 x l ) 


STF 


(3.2) 


The Newtonian and quadrupole-monopole contribution to the instantaneous loss under radiation reaction in the 
energy is : 


dE _ G ( 3 )ji, u (3 )ki 

— - ^5% { a N)I N 


dt 


( a n ) - ^- 5 -fjv^ ( a QM)4? W (ajv) • 


(3.3) 


The arguments of time derivatives of the momenta contain the contribution to the acceleration to be inserted in the 
respective terms. By inserting in the N ewto nian t erms the expressions (2.16) and (2.17) for v 2 and f 2 , respectively, 
and L(x) and Eqm{t , x) given by Eqs. ( 2.13 ) and ( 2.15 ) and the Newtonian expressions in the quadrupole-monopople 
terms, we obtain: 


dE 

dt 

dE\ 
dt ) 


dE 

UE\ 

dt 

\ dt ) N 

\ 

8 G 3 m 2 

)n 

15c 5 r 6 


2 G 3 m 4 

QM 

15c 5 Z 2 r 5 


dt J 


QM 


(2[xEr 2 + 2Gmfi 2 r+\lL 2 ) 


& ' { IN n cos (nx + Si) sin 2 m + 04(2 — 3 sin 2 Hi) 


(3.4) 

(3.5) 

(3.6) 


i=1 


< n—1 


The coefficients at are given by 

a\ = 3yiAr(—22Gmyi 2 r+ 17 L 2 ) , 


a 2 = 6 (-llG 2 m 2 /i 4 r 2 + 6 EL 2 /rr 2 + 5G?n/x 2 L 2 r- 51L 4 ) 
03 = — nAr(22Gmfi 2 r+ 51 L 2 ) , 

04 = 2 L 2 (-6E[ir 2 - 5Gm f j 2 r+ 39 L 2 ) . 


(3.7) 


By use of the true anomaly parametrization r(y), Eq. ( [2. If ) we find the energy loss in terms of \ alone. Then we 
pass to the complex parameter z = exp(*y) and we compute the averaged energy loss by use of the residue theorem, 
the only pole being at the origin: 
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dE\ _ /dE\ /dE\ 
dt / \ dt / N \ dt / 

dE\ = _ &m{-2mV> 
dt / N 15 c 5 L 7 


QM 


{U8E 2 L 4 + 732 G 2 m 2 n EL 2 + 425G 4 m> b ) , 


/dE\ GVm 3 (-2^) 3 / 2 ^ . 2 f ,^, 0 0 . 2 „ 

\-faj =- 30c 5 ^n -2^ Pism Ki cosSi + ° 2 ( 2 - 3 sin k*)] 

with the coefficients Ck given by: 

Ci = -iiA 2 (M8E 2 L 4 + 8936 G 2 m 2 p 3 EL 2 + 8335G 4 to 4 ^ 6 ) , 


C 2 = 708 E 3 L 6 + 10020 G 2 mV£ 2 T 4 + 18865G 4 toV £T 2 + 8316G b mV ■ 


^ 2^2 ,,3 7 t , 2 t 4 


' 4^,4 .6 Tp t 2 


^6 69 


(3.8) 

(3.9) 
(3.10) 


(3.11) 


B. Change in the magnitude of orbital angular momentum 


The instantaneous and averaged losses in L under the radiation reaction can be found in the same way as for the 
spin-spin case Q: 


dL ^ dJ ^ dSi ^ dS 2 

dt dt dt dt 


(3.12) 


The evolution under radiation reaction in the spin Si of the i th axisymmetric body (or approximately axisymmetric, 
with the deviation from axisymmetry of any post-Newtonian order), following |l3| was derived in |]24| by computing 
the integral of the moment of the reaction force (the sign swapped gradient of the Burke-Thorne potential) over the 
volume of each body. The radiation reaction will change the instantaneous orientation but not the magnitude of the 
spin vectors. However in [fl9| it was shown that there is no secular spin evolution under radiation reaction in the 2 PN 
order: 


dSj 

dt 


= 0 . 


(3.13) 


Den oting by ~ equalities modulo Burke-Thorne type terms, e.g. disregarding all terms which average out due to 
Eq.(|3.13) we find for the instantaneous loss in L under the radiation reaction: 


dL ^ dJ 
dt dt 

The leading order instantaneous loss of J under radiation reaction is: 

dJl _ 2G r m j(2).n T (3)kl 

dt 5 c 5 

The Newtonian and quadrupole-monopole terms decouple as: 

~AT = ~^ eyfe/ W b W)4f fc W) - ^ lN h \ a N)I [ N k \ a QM) + I ( N hl (a QM )lN )kl (^-n) 


(3.14) 


(3.15) 


(3.16) 


In the QM-terms we can use Cartesian coordinates ( x,y,z ) = r(cos ip, sin ip, 0) [equivalent to L = (0,0,1)]. Then the 
spins are expressed as S, =Sj(sin Ki cos xpi, sin m sin -i/’j, cos Ki). The computation yields the following instantaneous 
losses: 


t d 4 = (i d 4) + U d 4) 

dt \ dt) N \ dt J q M 
f dJ\ 8 G 2 mL , 2 -os 

L -dt) N =A^?r^ Er ’ 




2 G 


2 rm 3 


dt J QM 5c 5 p,L 3 r 7 


J2 Pi { b n cos (n\ + Si) sin 2 Ki + 64(2 — 3 sin 2 Ki) > , 


(3.17) 

(3.18) 

(3.19) 


s. n —1 
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with the coefficients bk given by 


b\ = 3^Ar(2Gm^i 3 Er 3 — fxEL 2 r 2 — 9Gm[i 2 L 2 r + 8 L 4 ) , 

b 2 = 3(2 G 2 m 2 n 5 Er 4 - 9 G 2 m 2 ^L 2 r 2 + 22 /iEL A r 2 + 10Gm/x 2 L 4 r - 19L 6 ) , 

63 = nAr(2Gmfj, 3 Er 3 + 3^EL 2 r 2 — 9Gmn 2 L 2 r — 24L 4 ) , 

64 = L 4 {—\8^Er 2 — 8Gm^x 2 r + 15 L 2 ) . 

A similar averaging procedure as in the case of the energy loss, gives here 


dL \ / dL\ / dL \ 

~dt ~ 


dt / n \ dt / q M 
/ dL\ AG 2 m(-2fj,E) 3 / 2 f2 , 1C ^ 2 3 , 

(*/„ = -- (14£i + 15G m " >' 

= G l| ' ''(-A' — y^p4DiSin 2 K,cos5i + D 2 (2-3sin 2 K j )] , 


where the coefficients Dk are given below 

Di = -6(62 E 2 L 4 + 211 G 2 m 2 v 3 EL 2 + 90G 4 mV) , 
D 2 = 252 E 2 L 4 + 1200 G 2 m 2 ii 3 EL 2 + 805G 4 mV . 


(3.20) 


(3.21) 

(3.22) 

(3.23) 


(3.24) 


C. Evolution of angles under radiation reaction 


cos 


The equations for the evolution under radiation reaction of the angles k, = cos 4 (Sj ■ L), (* = 1,2) and 7 = 

-it 


L (Si ■ S 2 ) were derived in 

d n 
dt 
d 

— COS Ki — I — COS K 


dt 


dt 


so 


d 

— COS K, 

dt 


ss 


d 

— COS K, 

dt 


QM 


(3.25) 

(3.26) 


The averaged SO terms were computed previously and given in |17j and the averaged SS terms in |]1J. To the 
expression derived in fly] we have added the required quadrupole-monopole term: 


— COS Ki 

dt 


QM 


~ i(Si - LcosKi) • 


QM 


(3.27) 


This is the new contribution to be computed in this subsection. For this purpose we rewrite the expression (3.27) as 
a function of the radial variables r(\) and \ : 


■=■ (Si — LCOS Ki) 

±j 


■ (A\ 


3 G 


2 -m3 


dt J q M 10 c 5 ^L 2 r 


+U 4 sin y sin 


• r, f ( 1 + 3j\ 

j sin Ki Pj sin 2k j < 2_ u n cos I n\ H-- —- 1 

7=1 n —1 ' 2 ' 


5i - Si 


■ U 5 COS 


Si - Sj 


)}• 


(3.28) 


Mi = —m 3 = ^,Ar(2^iEr 2 — 3 L 2 ) , 

M 2 = —2 L 2 (Gmfi 2 r + 3 L 2 ) , 

U 4 = 2{iAr(2/uEr 2 - 5 L 2 ) , 

m 5 = 2L 2 (AfiEr 2 + Gm/z 2 r - 5L 2 ) . (3.29) 

The residue theorem yields the averages: 
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(3.30) 


d'y 

dt 

diii 

dt 


/ dnj 

\ dt / so 


dKi 

dt 


dni 

dt 


SiS 2 


= 0 


dKi 

dt 


dni 

dt 


SS-self 


/ dKi ' 

\ dt / QM 


SO \ / SS-self 

given by Eq. (4.4) of |[7j 
given by Eq. (4.7) of |L9j 

given by Eq. (4.8) of ||T§ , 
3G 2 m 3 n(-2fiE) 3 / 2 A . o 


dt 


ss 


10 c 5 L 9 


j =i 


Vi cos 


where the coefficients V1-3 are: 



QM 


(h + <5j \ 

—) 


+ V 2 cos 



V) = 40 E 2 L 4 + 90 G 2 m 2 /i 3 EL 2 + 35 G 4 mV 
U 2 = 48 E 2 L 4 + U0G 2 m 2 fi 3 EL 2 + 70G 4 mV . 


(3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35) 

(3.36) 


For consisten cy in n otati on the replacements L —► L and Aq —► A should be carried out in Eq. (4.4) of |l^] . 

Equations ( ]3.30| )-(3.35) give the complete evolution under radiation reaction of the angles characterizing the relative 
orientation of the spin and orbital angular-momentum vectors up to second post-Newtonian order. 


IV. CONCLUDING REMARKS 


As the main result we have derived the second post-Newtonian order quadrupole-monopole contribution to the 
radiation reaction for the energy, the magnitude of the orbital angular momentum, and the angular variables charac¬ 
terizing the relative orientation of angular momenta of a coalescing binary system. This was possible by introducing 
a suitable angular average L together with the generalized true and eccentric anomaly parametrizations for the 
quadrupole-monopole perturbati on of the Keplerian motion. 

The secular losses of E, Eqs. (3.8)-( 3.10 ) and of L, Eqs. (|3.21 )-( T23| ), complementing the corresponding spin-orbit 
terms of Ref. [fl7| (with the replacements L —> L and Ao —* A), the spin-spin terms of Ref. [Q, and the post-Newtonian 
terms of Ref. |pd| give the total radiation reaction up to second post-Newtonian order. 

All contributions to the secular evolutions of the angles Kj and 7 were listed in the Sec III.C. They are of either 
spin or quadrupolar origin. 

Finally, we compare our result for the energy loss w ith th e one given by Poisson in Ref. [ p0| for circular orbits. For 
this we need to take the circular orbit limit of Eqs. ( |3.8| )-(3.11). By the circular orbit limit in Ref. [|o| is meant a 
perturbation of a circular Keplerian orbit. This, however, cannot be immediately imposed in our formalism, relying 
on heavy use of conserved quantities of the perturbed motion. It is certainly impossible to impose the circularity 
conditions for the unperturbed orbit as En = E — Eqm = —Gmy/2ro and L 2 N = L 2 (y) = (L + SL) 2 = Gmfi 2 ro 
(ro being the radius of the unperturbed orbit) due to the fact that Eqm and 5L are not constants, while E , L and 
ro are. This is related to the remark of ( 2 ^] that the orbits cannot be circular in the strict sense. However, we can 
impose the circularity conditions in an average sense. As the angular average of 5L vanishes, we obtain: 


En = E — Eqm = — 


Gmyi 


2rn 


L 2 n = L 2 = Gmp ro . 


(4.1) 


These values inserted in the angular average of Eqm (which is already of second PN order) give 


Gm 3 p \ ^ / 

e QM = ——2^ Pi i 3 

i =1 


Ar 3 


cos /y 


- 1 ) 


(4.2) 


Expressing E and L from Eqs. © and 
we find the radiative energy loss for the a 


2 ) and inserting them into the expressions 
rove defined circular orbits: 


(4 


-(3.11) for the energy loss 
























32 G 4 my 
5c 5 rQ 


TH? ^ 

1 - 6— ^ (3 cos 2 Kt — l) 

r o i=l 


(4.3) 


Employing the relation between the radius of the unperturbed circular orbit ro and the average velocity of the 
perturbed orbit (v), deducible from the expressions given in Ref. |2(j| : 


r o = 


Gm 

w 


1 “ ( 3c0s2 K i - !) 


we find that the averaged energy loss for circular orbits agrees with Eq. (22) of Ref. 


(4.4) 
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